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ABSTRACT 

We analyze a set of volume limited samples from SDSS DR12 to quantify the 
degree of inhomogeneity at different length scales using Shannon entropy. We find 
that the galaxy distributions exhibit a higher degree of inhomogeneity as compared 
to a Poisson point process at all length scales. Our analysis indicates that signa¬ 
tures of inhomogeneities in the galaxy distributions persist at least upto a length 
scale of 120/1“'Mpc. The galaxy distributions appear to be homogeneous on a scale 
of 140/1“* Mpc and beyond. Analyzing a set of mock galaxy samples from a semi ana¬ 
lytic galaxy catalogue from the Millennium simulation we find a scale of transition to 
homogeneity at ~ 100/i“*Mpc. 

Key words: methods: numerical - galaxies: statistics - cosmology: theory - large scale 
structure of the Universe. 


1 INTRODUCTION 

Homogeneity and isotropy on sufficiently large scales is an 
assumption which is fundamental to our current understand¬ 
ing of the Universe. This assumption which is popularly 
known as the Cosmological principle implies that the prop¬ 
erties of the Universe is same for all observers irrespective 
of their locations and the directions at which they are look¬ 
ing at. The assumption has also an aesthetic appeal of its 
own as it indirectly asserts us that the same physical laws 
must apply throughout the Universe given it is homogeneous 
and isotropic to start with. One can not prove the Cosmo¬ 
logical principle in a strictly mathematical sense but fortu¬ 
nately one can test it against the presently available cosmo¬ 
logical observations. The existence of the cosmic microwave 
background radiation and its near uniform temperature over 
the entire sky provides by far the best conclusive evidence 
in favour of isotropy (Penzias & Wilson 1965; Smoot et al. 
1992; Fixsen et al. 1996). Other important supporting evi¬ 
dences come from the isotropy in angular distributions of 
radio sources (Wilson & Penzias 1967; Blake & Wall 2002), 
isotropy in the X-ray background (Peebles 1993; Wu et al. 
1999; Scharf et al. 2000) and the isotropy in the distribu¬ 
tion of neutral hydrogen as shown by lyman-a transmit¬ 
ted flux (Hazra & Shafieloo 2015). But mere existence of 
isotropy around us does not automatically guarantees ho- 
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mogeneity of the Universe. Spatial homogeneity can be as¬ 
serted from isotropy only when the later is confirmed around 
each point with the hypothesis that the matter distribu¬ 
tion is a smooth function of position (Straumann 1974; 
Sylos Labini & Baryshev 2010). 

We see structures in the present Universe starting from 
planets, stars and galaxies to groups, clusters and super¬ 
clusters spanning a wide range of length scales. The present 
Universe looks highly inhomogeneous on small scales and 
it is important to test if such inhomogeneities continue to 
persist on large scales. Large scale inhomogeneities have sev¬ 
eral important implications for cosmology and could pose a 
serious challenge to the standard cosmological framework. 
Inhomogeneities, through the backreaction mechanism can 
provide an alternate explanation of a global cosmic acceler¬ 
ation without requiring any additional dark energy compo¬ 
nent (Buchert & Ehlers 1997; Buchert 2001; Schwarz 2002; 
Kolb et al. 2006; Paranjape 2009; Kolb et al. 2010; Ellis 
2011). A large void can also mimic an apparent acceleration 
of expansion (Tomita 2001). Fortunately such models can be 
constrained with observations such as SNe, CMB, BAO and 
measurements of angular diameter distance and Hubble pa¬ 
rameter (Zibin et al. 2008; Clifton et al. 2008; Biswas et al. 
2010; Clarkson 2012; Larena et al. 2009; Valkenburg et al. 
2014). However if our Universe is inhomogeneous on large 
scales, the currently available inhomogeneous cosmological 
models suggest that the effects of inhomogeneities on ob¬ 
served quantities are non-negligible and more precise obser- 


© 2015 The Authors 


2 Pandey, B. and Sarkar, S. 


vations will not be properly analyzed unless inhomogeneities 
are taken into account (Bolejko et al. 2011). 

Galaxy surveys map the three dimensional distribution 
of galaxies in the Universe. If our Universe is homogeneous 
then the statistical properties of galaxy distributions in a 
finite volume should be independent of the location of that 
volume in the Universe. One simple tool to characterize the 
statistical properties of galaxy distributions is the two point 
correlation function (Peebles 1980) which measures the ex¬ 
cess probability of finding a pair of points separated by a 
distance r compared with a random Poisson process. The 
two point correlation function on small scales 0.1/r“*Mpc < 
r< lOh^'Mpc, is well described by a power law of the form 
^(r) = with correlation length ro ~ 5-6/2^*Mpc and 

slope 7 ~- 1.7- 1.8 (Hawkins et al. 2003; Zehavi et al. 2005). 
^(r) vanishes at scales > 20/P*Mpc which is consistent with 
large scale homogeneity. However the problem with corre¬ 
lation function analysis is that it assumes a mean density 
on the scale of survey which is not a defined quantity be¬ 
low the scale of homogeneity. Most of the statistical tests 
of homogeneity carried out so far are based on the count 
in spheres method where the number counts n(< r) within 
a sphere of radius r and its scaling with r is used to test 
the transition scale to homogeneity. If the scaling expo¬ 
nent approaches a value ~ 3 on some scale it marks the 
transition to homogeneity indicating that the distribution 
is homogeneous on and above that scale (Hogg et al. 2005; 
Sylos Labini 2011a; Scrimgeour et al. 2012). Fractal anal¬ 
ysis (Martinez & Jones 1990; Coleman & Pietronero 1992; 
Borgani 1995; Bharadwaj et al. 1999; Yadav et al. 2010) 
uses the scaling of different moments of n(< r) to char¬ 
acterize the scale of homogeneity. Various studies carried 
out with these methods claim to have found a transition 
to homogeneity on sufficiently large scales 70- 150/i“*Mpc 
(Martinez & Coles 1994; Guzzo 1997; Martinez et al. 1998; 
Bharadwaj et al. 1999; Pan & Coles 2000; Kurokawa et al. 
2001; Hogg et al. 2005; Yadav et al. 2005; Sarkar et al. 
2009; Scrimgeour et al. 2012) whereas there are stud¬ 
ies which claim the absence of any such transition 
out to the scale of the survey (Coleman & Pietronero 
1992; Amendola & Palladino 1999; Sylos Labini et al. 2007, 
2009a,b; Sylos Labini 2011b). The present generation galaxy 
surveys (SDSS, York et al. 2000; 2dFGRS, Colles et al. 
2001) have now mapped the distribution of millions of galax¬ 
ies across billions of light years providing an unique oppor¬ 
tunity to test the assumption of homogeneity on large scales. 
The Sloan Digital Sky Survey (SDSS) has reached its final 
stage (Data Release 12) (Alam et al. 2015) now encompass¬ 
ing more than one-third of the entire celestial sphere. This 
provides us with galaxy distributions in enormous volumes 
of the Universe for which the statistical properties of galaxy 
distributions can be analyzed and homogeneity of the Uni¬ 
verse can be tested with greater confidence. Pandey (2013) 
introduce a method based on Shannon entropy (Shannon 
1948) for characterizing inhomogeneities and applied the 
method on some Monte Carlo simulations of inhomogeneous 
distributions and N-body simulations which show that the 
proposed method has great potential for testing the large 
scale homogeneity in galaxy redshift surveys. 

In the present work we employ the method proposed 
by Pandey (2013) to analyze a set of volume limited sam¬ 
ples from the Sloan Digital Sky Survey Data Release Twelve 


(SDSS DR12). We quantify the inhomogeneities in the 
galaxy distribution using Shannon entropy and test if there 
is a scale of transition to homogeneity. We also analyze a 
set of mock galaxy samples from a semi analytic galaxy cat¬ 
alogue (Guo et al. 2011) derived from the Millennium Run 
simulation (Springel et al. 2005) to compare theoretical pre¬ 
dictions with observations. 

A brief outline of the paper follows. We briefly describe 
our method in Section 2, describe the Monte Carlo simula¬ 
tion, SDSS and Millennium data in Section 3 and present 
the results and conclusions in Section 4. 

We have used a ACDM cosmological model with = 
0.3, Dao = 0-7 and h = \ throughout. 


2 METHOD OF ANALYSIS 

Pandey (2013) propose a method based on the Shannon 
entropy to study inhomogeneities in a 3D distribution of 
points. Shannon entropy (Shannon 1948) is originally pro¬ 
posed by Claude Shannon to quantify the information con¬ 
tent in strings of text. It gives a measure of the amount 
of information required to describe a random variable. The 
Shannon entropy for a discrete random variable X with n 
outcomes [xj : i = is a measure of uncertainty denoted 

by H(X) defined as. 


n 

H{X) = - ^ p{xi) log p{xi) (1) 

where p{x) is the probability distribution of the random 
variable X. 

Given a set of N points distributed in 3D we consider 
each of the points as center and determine ra,(< D the 
number of other points within a sphere of radius r as. 


N 

nA< D = 20U-|x/-X;I) 
Al 


( 2 ) 


where 0 is the Heaviside step function and Xj and xj are the 
radius vector of and p points respectively. To avoid any 
edge effects we discard all the centers which lie within a dis¬ 
tance r from the survey boundary. The number of available 
centers will decrease with increasing r for any finite volume 
sample. We define a separate random variable W for each 
radius r which has M{r) possible outcomes each given by, 
fir= — with the constraint yPP fir = L The Shan- 

non entropy associated with the random variable X,. can be 
written as, 


Hr 


M(r) 

- 2 fur log fi,r 
i=\ 

M{r) 

log)^ «/(< D)- 

i=\ 


a(< r) log(n,(< r)) 
Zrj[^n,«r) 


( 3 ) 


Where the base of the logarithm is arbitrary and we choose it 
to be 10. fi^r will have the same value for all the centers 
when ra;(< r) is same for all of them. This is an ideal situation 
when the spheres of radius r around each of the M(r) cen¬ 
tres contain exactly the same number of points. This max¬ 
imizes the Shannon entropy to {Hr)max = log M(r) for radius 
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r. We define the relative Shannon entropy as the ratio of the 
entropy of a random variable to the maximum possible 
entropy {Hr)max associated with it. The relative Shannon en¬ 
tropy at any r quantifies the degree of uncertainty in 

the knowledge of the random variable Xy. The distribution 
of fr become completely uniform when = 1 is reached. 

In the present scenario the random variables Xy are not 
mutually independent as the measurements around different 
centers at each radius r are carried out on the same finite 
volume sample and multiple spheres around the centers can 
overlap. The random variables may have extra correlations 
if the points are clustered or distributed in an intrinsically 
inhomogeneous way. These extra correlations increase the 
mutual information stored in the random variables Xy. The 
decrease of information in Xy with increasing r is evident 
irrespective of a homogeneous/inhomogeneous distribution 
but it would diminish differently depending on the nature 
and the degree of inhomogeneity present in the distribution. 
When inhomogeneities are present, the relative Shannon en¬ 
tropy would show a departure from the values 1. Larger de¬ 
parture clearly indicates greater degree of inhomogeneity. 
When the Shannon entropy Hy attains its maximum value 
(Hy)max corresponding to that radius r this ratio levels up 
with 1. This indicates absence of inhomogeneity beyond that 
scale and may be considered as the scale of transition to 
homogeneity. But this should be taken with some caution 
as overlap between spheres around different centres forcibly 
brings down the inhomogeneities by increasing the relative 
Shannon entropy at each length scales. Overlapping become 
progressively more important at larger radii. The number 
of available centres at a radius gradually decreases with in¬ 
creasing r and the centres preferentially migrate towards the 
centre of the volume analyzed. The available centres are fi¬ 
nally confined to a smaller region when the largest spheres 
are considered. This confinement bias (Pandey 2013) may fi¬ 
nally force the relative Shannon entropy to 1 at some radius 
r when the spheres completely overlap with each other. This 
corresponds to a situation where the spheres around differ¬ 
ent centres show a near uniform number count mimicking 
a real transition to homogeneity. The degree of overlap as 
well as the confinement bias would depend on the shape of 
the volume and nature of the distribution to some extent. A 
large volume ensures spheres upto a larger length scales and 
hence may help us to detect the scale of homogeneity if the 
transition occurs before the scale where the confinement bias 
completely dominates the statistics. For the present analysis 
one may also consider non-overlapping spheres of different 
radii. But the statistics becomes too noisy due to very small 
number of independent spheres at progressively larger radii 
which consequently prohibits us to address the issue of ho¬ 
mogeneity on large scales. 

The observed inhomogeneities are the combined out¬ 
come of gravitational clustering, Poisson noise and any in¬ 
trinsic inhomogeneity present in the distribution. Inhomo¬ 
geneity due to Poisson noise rapidly decreases with r. The 
inhomogeneity due to clustering would also decrease with r 
albeit differently depending on the strength and nature of 
clustering. Overlap causes a further suppression in the in¬ 
homogeneities at all scales. A complete discussion on the 
various sources of inhomogeneity and the factors affecting 
the value of the relative Shannon entropy can be found in 
Pandey (2013). 


3 DATA 

3.1 MONTE CARLO SIMULATIONS 

Application of the method discussed in Section 2. involves 
analysis of galaxy distributions in finite regions of the Uni¬ 
verse. Construction of volume limited samples from galaxy 
surveys over different redshift and magnitude ranges would 
result in samples having different volumes and numbers of 
galaxies. One may also prefer to analyze volume limited 
galaxy samples having different geometries. Keeping these 
possibilities in mind we would like to first investigate how 
the parameters like number density, volume and geometry 
affect the relative Shannon entropy of a 3D spatial distribu¬ 
tion of a point process. 

To study the possible systematic effects we generate dif¬ 
ferent sets of Monte Carlo realizations of the homogeneous 
Poisson point process. We consider the following cases: (i) 
number density is allowed to vary keeping the volume and 
geometry fixed, (ii) volume is allowed to vary keeping the 
number density and geometry fixed and (iii) geometry is 
allowed to vary keeping both the number density and vol¬ 
ume fixed. We simulate 10 Monte Carlo realizations for each 
density, volume and geometry considered. Table 1, Table 2 
and Table 3 summarize the properties of the Monte Carlo 
simulations generated to test the effects of number density, 
volume and geometry respectively. 


3.2 SDSS DR12 DATA 

The Sloan Digital Sky Survey (SDSS) which started its op¬ 
erations with a 2.5m telescope (Gunn et al. 1998, 2006) in 
2000 has now reached its final stage after carrying out imag¬ 
ing and spectroscopy over one third of the Celestial sphere. 
SDSS DR12, the final data release of SDSS (Alam et al. 
2015) contains all data taken by all phases of the SDSS 
through 14th July, 2014. DR12 incorporate some significant 
improvements to the spectrophotometric flux calibration. It 
now contain optical spectroscopy of 2401952 galaxies and 
477161 quasars representing a ~40% increase in Baryon Os¬ 
cillation Spectroscopic Survey (BOSS) spectra over the pre¬ 
vious data release DRIO (Ahn et al. 2014). 

Our present analysis is based on SDSS DR12 data. 
We have used the Main Galaxy Sample which comprises of 
galaxies brighter than a limiting r-band Petrosian magnitude 
17.77. The target selection algorithm of the Main Galaxy 
Sample is detailed in Strauss et al. (2002). We downloaded 
the data from the Catalog Archive Server (CAS) of SDSS 
DR12 using a Structured Query Language (SQL) search. 
We identify a contiguous region spanning 135 < a < 225 and 
0 < h < 60 where a and 6 are the equatorial co-ordinates. We 
construct a set of volume limited samples by restricting the 
extinction corrected r-band Petrosian apparent magnitude 
to the range 14.5 < my < 17.5. The r-band absolute magni¬ 
tude ranges, the corresponding redshift ranges, number of 
galaxies, total volume, number density, mean inter-particle 
separation and the size of the largest sphere that can com¬ 
pletely fit inside each of the volume limited samples are listed 
in Table 4. 


MNRAS 000, 1-12 (2015) 




4 Pandey, B. and Sarkar, S. 


Table 1. This shows the properties of the Monte Carlo simulations of homogeneous Poisson process generated inside a cube maintaining 
exactly the same volume while varying the number density. The size of the largest sphere tabulated here are the averages of the same 
over 10 different realizations used in each case. 


Sample 

name 

Box Length 
{h~^ Mpc) 

Volume 
{h~^ Mpc)^ 

No. Density 

(/r‘ Mpc)-3 

Number of 
points 

Largest Sphere 
{h~^ Mpc) 

density 1 

200 

8x10^ 

5.0x10-’* 

4000 

95 

density 2 

200 

8x10^ 

1.0x10-3 

8000 

96 

density 3 

200 

8x10^ 

2.0x10-3 

16000 

96 

density 4 

200 

8x10'’ 

4.0x10-3 

32000 

97 


Table 2. Same as Table 1 but here the volumes of the samples are allowed to change while keeping the number density and geometry 
same. 


Sample Box Length Volume No. Density Number of Largest Sphere 
name (/ 2 “^Mpc) (/T^Mpc)^ (/i“^Mpc)“^ points (/ 2 “^ Mpc) 


volume 1 

100 

LOxlO* 

5.0x10-3 

4999 

48 

volume 2 

125 

1.95x10'’ 

5.0x10-3 

9765 

60 

volume 3 

150 

3.38x10** 

5.0x10-3 

16874 

72 

volume 4 

175 

5.36x10*’ 

5.0x10-3 

26796 

85 


Table 3. Same as Table 1 but here the geometry of the samples are allowed to change while keeping the number density and volume 
same. 


Geometry 

Boundary 

Volume 
(ft-* Mpc)3 

No. Density 
(/i“‘ Mpc)“^ 

Number of 
points 

Largest Sphere 
{h~^ Mpc) 

Survey like 

0.02 < j < 0.08 

135 < O'<225 
0<6<60 

5.82x10** 

2.0x10-3 

11633 

78 

Sphere 

radius:111.7 Mpc 

5.82x10** 

2.0x10-3 

11633 

107 

Cube 

side:179.8 ft-* Mpc 

5.82x10** 

2.0x10-3 

11633 

87 


Table 4. This summarizes the properties of our volume limited samples from SDSS. 


Galaxy 

Sample 

Absolute 
magnitude range 

Redshift 

range 

Number of 

Galaxies 

(N) 

Volume of 
the region 
(ft-* Mpc)3 

Number 

density 

Mpc)“^ 

Mean 
separation 
(ft-* Mpc) 

Largest 

Sphere 
{h~^ Mpc) 

sample 1 

-20.5 >Mr> -22 

0.067 <z< 0.125 

51804 

1.86x10** 

2.77x10-3 

7.1 

83 

sample 2 

-20.75 >Mr> -22 

0.067 <z< 0.138 

48992 

2.61 xlO’ 

1.87x10-3 

8.1 

103 

sample 3 

-21 >Mr> -22 

0.067 <z< 0.153 

40263 

3.58 xl03 

1.12x10-3 

9.6 

122 

sample 4 

-21.25 > Mr >-22 

0.067 <z< 0.168 

30196 

4.84x10** 

6.23x10-’* 

11.7 

144 

sample 5 

-21.5 > Mr >-22 

0.067 <z< 0.185 

19063 

6.45x10** 

2.95x10-’* 

15 

167 

sample 6 

-21.7 > Mr >-22 

0.067 <z< 0.199 

10567 

8.05 xlO’ 

1.31x10-’* 

19.6 

185 
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Red shift (z) 


Figure 1. This shows the definition of our six volume limited samples (Table 4) from SDSS in redshift-absolute magnitude plane. The 6 
overlapping boxes delineate 6 samples used in this analysis. The two curves define the boundary corresponding to the apparent magnitude 
limit imposed. 


3.3 MILLENNIUM DATA 

Semi analytic models (White & Frenk 

1991; Kauffmann, White & Guiderdoni 1993; 

Kauffmann & White 1993; Kauffmann 1996; Cole et al. 
1994, 2000; Somerville & Primack 1999; Baugh et al. 1998; 
Benson et al. 2002; Springel et al. 2005; Guo et al. 2011) 
provide a very powerful tool to study galaxy formation and 
evolution. Galaxy formation and evolution involve many 
physical processes such as gas cooling, star formation, 
supernovae feedback, metal enrichment, merging and mor¬ 
phological evolution. The semi analytic models parametrise 
the physics involved in terms of simple models following 
the dark matter merger trees over time. The models 
provide the statistical predictions of galaxy properties 
at some epoch and the precision of these predictions are 
directly related to the accuracy of the input physics. In the 
present analysis we use a semi-analytic galaxy catalogue 
generated by Guo et al. (2011) from the Millennium Run 
simulation (Springel et al. 2005) who updated the previ¬ 
ously available galaxy formation models (Springel et al. 
2005; Croton et al. 2006; De Lucia & Blaizot 2007) with 
improved versions. The spectra and magnitude of the 
model galaxies were computed using population synthesis 
models of Bruzual & Chariot (2003). We map the galaxies 
to redshift space using their peculiar velocities and then 
identify regions which have the same geometry as our SDSS 
samples. We then apply the same magnitude cuts as those 
used for the actual data and extracted the same number of 
galaxies as in our SDSS samples. For each of the volume 


limited samples used in our analysis we generate three mock 
galaxy samples from the semi analytic galaxy catalogue. 


4 RESULTS AND CONCLUSIONS 

In the top left Figure 2 we show the variations of the rel¬ 
ative Shannon entropy jjpj — with distance r for some ho¬ 
mogeneous Poisson point processes generated inside a cu¬ 
bic box with different number densities (Table 1). The rela¬ 
tive Shannon entropy versus r curves show a characteristic 
behaviour in all cases. At the smallest radius the relative 
Shannon entropy starts from a value which is very close to 
1. It suddenly drops at a particular length scale and then 
starts increasing again as the length scale increases. The 
behaviour of Shannon entropy at the smallest radius and 
the presence of a characteristic dip in it on a certain length 
scale can be understood as follows. Any distribution having 
a finite number density will look nearly uniform below the 
scale of mean inter-particle separation corresponding to that 
distribution. As this characteristic length scale is surpassed 
the irregularities or any deviations from homogeneity can 
be clearly seen. The dips in the set of curves in all panels 
of Figure 2 are located at the mean inter-particle separa¬ 
tions of the respective samples. While testing the effects of 
number density we have kept unaltered all the other fac¬ 
tors such as volume and the geometry of the samples. As 
the density of the sample drops the characteristic dip in 
the relative Shannon entropy shifts towards larger length 
scales (top left panel of Figure 2) simply resulting from an 
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Figure 2. The top left, top right and lower middle panel show the systematic effects of number density, volume and geometry of samples 
respectively on the relative Shannon entropy of a Poisson process at different length scales. The error bars shown here are the 1 — cr 
variations from the 10 realizations used in each case. 


increase in the mean inter-particle separation. The differ¬ 
ences in the amplitudes of the relative Shannon entropy for 
Poisson samples with different densities result from the dis¬ 
creteness noise or Shot noise due to Poisson fluctuations. As 
the density of the sample decreases the Poisson noise shoots 
up at each length scale causing a decrease in the relative 
Shannon entropy and increase in the degree of inhomogene¬ 
ity at the corresponding scales. The points are uncorrelated 
in case of a Poisson point process. The Void Probability 
function (White 1979) for a homogeneous Poisson distribu¬ 
tion is where A is the intensity of the Poisson process 
and V is the volume. It is expected that such a distribution 
would become homogeneous on a scale r ~ A^i, which is a 
measure of the average size of voids in the distribution. This 
is closely related to the mean inter-particle separation and 
one would expect a Poisson distribution to become homo¬ 
geneous above this length scale. But one can see in the top 
left panel of Figure 2 that the relative Shannon entropy of 
all the Poisson point processes continue to show a departure 
from 1 upto some distance even after this length scale. This 
is due to the Poisson noise which rapidly decreases with 
increasing r. Further the spheres around different centres 


overlap with each other which reduces the inhomogeneity 
increasing the relative Shannon entropy. Overlapping also 
becomes progressively important at larger length scales in 
presence of confinement bias. The Poisson noise and over¬ 
lapping affect the relative Shannon entropy in opposite man¬ 
ner in which the former enhances whereas the later reduces 
inhomogeneities at each length scales. It may be noted that 
the Poisson noise dominates at smaller scales whereas over¬ 
lapping dominates at larger scales. This eventually allows 
the overlapping to wipe out the inhomogeneities introduced 
by the Poisson noise. It can be seen in the top left panel of 
Figure 2 that the relative Shannon entropy eventually level 
up with 1 for all the Poisson point processes but at different 
length scales. The lowest density sample having the largest 
Poisson noise at each length scale become homogeneous on 
the largest length scale among the samples analyzed. We 
have used here an extreme variation in density for this test, 
the highest density samples being 8 times denser than the 
lowest density samples (see Table 1). We also test the affects 
of the size of volume occupied by the distribution. We dras¬ 
tically vary the volume of the samples keeping their number 
density and geometry unchanged (see Table 2). The results 
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Figure 3. Different panels shows the relative Shannon entropy at different length scales for the first four among six SDSS samples 
described in Table 4. The sample names and the corresponding redshift and absolute magnitude ranges are displayed at the top right 
corner of each panel. Each panel also show together the results for the corresponding mock galaxy samples and Poisson samples for 
comparison. The results for sample 5 and sample 6 are shown separately in Figure 4 as mock samples for them could not be prepared 
from the Millennium simulation due to size restrictions. The error bars shown here are the 1 — cr variations from the 10 and 3 mock 
samples for the Poisson process and Millennium simulation respectively. The error bars are very tiny and nearly invisible here. The 
relative Shannon entropy for all the samples are computed at intervals of Ift^'Mpc (as in Figure 2) but for clarity we show here the 
results only at intervals of Mpc. 


are shown in the top right panel of Fignre 2. We can see 
that the relative Shannon entropy are largely unaffected by 
the change in volume. It shows a very small increase in rela¬ 
tive Shannon entropy nearly on all scales, the changes being 
noticeable only when the volume is increased by a factor of 
~ 6. This small change in Shannon entropy can be attributed 
to the cosmic variance as the number of centers available at 
a particular radius increases with the volume. The largest 
length scale that can be probed with these samples are dif¬ 
ferent as they cover different volumes (Table 2). It is in¬ 
teresting to note that all the samples become homogeneous 
at nearly the same length scales 30Mpc despite the fact 
that the samples with different volumes probe upto different 
length scales and have different number of centers at each 
radii. Although these samples cover different volumes they 


have exactly the same number density ensuring the same 
contribution from Poisson noise at each length scales. This 
clearly demonstrates that the Poisson noise or discreteness 
noise governed by the number density of the samples plays 
an important role in deciding the degree of inhomogeneity 
and transition to homogeneity. Finally we test the effects of 
geometry of the volumes on the relative Shannon entropy. 
We use a survey like region, a cubic box and a spherical re¬ 
gion all having the same number density and volume (see 
Table 3). The results are shown in the lower middle panel 
of Figure 2. We see that the degree of inhomogeneity and 
the transition scale to homogeneity both are unaffected by 
the change in geometry of the samples. The results of these 
tests clearly indicate that the relative Shannon entropy of a 
distribution is most sensitive to the number density of the 
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Figure 4. This shows the relative Shannon entropy at different length scales for sample 4 together with the same for two larger samples 
namely sample 5 and sample 6 described in Table 4. Despite having different volumes and probing upto different scales the results from 
all three samples merge together at ~ 140/i“* Mpc. The offset between the curves are primarily due to differences in number densities and 
clustering which cease to exit beyond 140/2“* Mpc . As in Figure 3 we also show here the results at intervals of 5/2“* Mpc for clarity. 
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Figure 5. This shows the number of available centres at different length scales for each of the SDSS samples analyzed. The length scales 
where there are only 100 centres for each of the SDSS samples are shown with arrows at the bottom. We show here the results at intervals 
of 1/r* Mpc. 
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Figure 6. This shows the relative Shannon entropy at different length scales for sample 4 together with the same for sample 5 and 
sample 6 after appropriate cuts are applied to them to match their sizes with sample 4. 



Figure 7. This shows the number of available centres at different length scales for sample 4 together with the same for sample 5 and 
sample 6 after appropriate cuts are applied to them to match their sizes with sample 4. 
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distribution and nearly insensitive to the volume and the 
geometry of the samples. 

In Figure 3 we show the relative Shannon entropy at dif¬ 
ferent length scales for the first 4 SDSS samples described 
in Table 4 together with that from their mock counterparts 
from the Millennium simulation and Poisson point processes. 
These four different SDSS volume limited samples are con¬ 
structed by keeping the brighter magnitude limit fixed at 
-22 while gradually shifting the fainter magnitude limit from 
-20.5 to -21.25 in steps of -0.25. This ensures minimal dif¬ 
ferences in the clustering properties of the resulting samples 
while allowing us to probe different length scales due to their 
different volume coverages. The mock samples drawn from 
the Millennium simulation and Poisson processes have iden¬ 
tical number density, volume and geometry as their respec¬ 
tive actual SDSS counterparts. We see a characteristic dip 
in the relative Shannon entropy curve for all the samples at 
a length scale which corresponds to their mean inter-particle 
separations (Table 4). Further in all the panels of Figure 3 
we see that the galaxy samples from both the SDSS and 
the Millennium have a higher degree of inhomogeneity as 
compared to their mock Poisson samples. These differences 
indicate that the galaxy samples must have extra sources of 
inhomogeneity other than the Poisson noise resulting purely 
from the discrete nature of the distributions. It is well known 
that galaxy distributions exhibit clustering which can be 
quantified by their correlation functions. Clustering acts as 
an extra source of inhomogeneity for the galaxy samples. 
Besides there could be an intrinsic inhomogeneity built in 
the distribution. The inhomogeneity from clustering and in¬ 
trinsic inhomogeneous nature if any combine together with 
that from Poisson noise to produce the total inhomogeneities 
in the distributions. Irrespective of their nature inhomo¬ 
geneities are expected to decrease with increasing length 
scales in presence of overlapping and confinement bias. Con¬ 
sequently these together makes it difficult to distinguish a 
real transition to homogeneity from an induced one. The 
transition can be confirmed easily if it takes place at a length 
scale where confinement bias does not dominate the statis¬ 
tics. For example in Figure 3 the mock Poisson samples cor¬ 
responding to four SDSS samples show a transition to ho¬ 
mogeneity in the range 30 - 50/;“* Mpc. Although the tran¬ 
sitions are affected by Poisson noise on small scales still we 
can make a decision on the existence of a genuine transition 
to homogeneity. Galaxy distributions are far from a random 
Poisson distribution and if a scale of transition to homo¬ 
geneity exist there it would possibly occur on sufficiently 
large scales where the confinement bias may overpower the 
inhomogeneities. Interestingly one can reduce the impact of 
confinement bias by simply choosing a volume which are less 
symmetric than a sphere or a cube. Although we have men¬ 
tioned earlier that the geometry of the analyzed volume does 
not affect the relative Shannon entropy as shown in Figure 2, 
it could actually affect the statistics if the distribution is in¬ 
homogeneous nearly upto the largest length scale that can 
be probed with that volume. A volume like a rectangular 
slab whose two dimensions are significantly larger than the 
third one would make it sure that the available centres are 
spread across nearly the entire volume and would therefore 
never allows a situation where the spheres completely over¬ 
lap with each other even at the largest radius. Thus if the 
number counts around different centres differ even slightly 


the relative Shannon entropy will show a deviation from 1. 
The geometry of our SDSS samples are not exactly slab like 
but definitely has a lesser symmetry than a spherical or cubic 
volume. This may enable us to detect the signature of any 
inhomogeneities present at the largest length scales probed 
by these samples. Although overlapping would be present at 
each scale, it may not be able to erase the inhomogeneities 
completely even at the largest scale due to the lesser im¬ 
pact of the confinement bias. Any residual inhomogeneity 
detected at the largest length scales are less likely to borne 
out of pure Poisson noise. At the largest length scale the 
numbers of spheres could be drastically small so we decided 
to consider only upto the length scales where this number 
reduces to 100. In Figure 5 we show how the number of cen¬ 
ters vary with length scales for all the six SDSS samples 
analyzed here. We can see that different samples probe upto 
different length scales and the scale where the number of 
available centres reduces to 100 for sample 1, sample 2, sam¬ 
ple 3 and sample 4 are 83/;“*Mpc, 102/;“* Mpc, 119/;“* Mpc 
and 138/;“* Mpc respectively. In the top left panel of Figure 3 
the relative Shannon entropy for the SDSS sample 1 clearly 
shows a deviation from unity at 80/;“* Mpc signaling some 
residual inhomogeneities on that scale. The results of the 
mock Poisson samples shown in Figure 3 indicate that for 
the galaxy samples having number densities comparable to 
our SDSS samples the contribution from Poisson noise can 
be safely neglected beyond 50/;“* Mpc. Further sample 1 be¬ 
ing the densest of the SDSS samples analyzed here would be 
least affected by the Poisson noise at all length scales. This 
indicates that the galaxy distribution is not homogeneous on 
80/;“* Mpc. The results from SDSS sample 2, sample 3 and 
sample 4 are shown in the top right, bottom left and bottom 
right panel of Figure 3 respectively, sample 2, sample 3 and 
sample 4 are 1.4, 2.4 and 4.4 times less denser than sample 
1 (Table 4) and are expected to have a greater contribution 
from Poisson noise at all length scales. Interestingly the rel¬ 
ative Shannon entropy for sample 1, sample 2, sample 3 and 
sample 4 are nearly same at 80/;“* Mpc. As sample 2, sample 
3 and sample 4 probe progressively larger length scales we 
see that the relative Shannon entropy increases further with 
increasing length scales. It may appear as if the reduction in 
inhomogeneity with increasing length scales as revealed by 
progressively larger samples could result from a greater de¬ 
gree of overlap between the spheres. But the samples probing 
larger scales have larger volumes and for each of them we 
restrict our attention only upto a scale where the number of 
centers reduces to 100. So the degree of overlap would not 
vary across the samples given we limit our attention upto 
this specihc length scale. All the four panels of the Figure 3 
show that the relative Shannon entropy change very slowly 
on large scales indicating that the degree of overlap do not 
change rapidly for our sample geometries possibly due to 
a less aggressive confinement bias. It may be noted further 
that the values of relative Shannon entropy for sample 2, 
sample 3 and sample 4 are nearly same at 102/;“* Mpc and 
they are also same at 119/;“* Mpc for sample 3 and sample 
4. These together clearly demonstrate the decreasing impor¬ 
tance of Poisson noise with increasing length scales and indi¬ 
cate that any residual inhomogeneities at the largest length 
scale could be safely attributed to the presence of genuine in¬ 
homogeneity in the distribution at the corresponding length 
scale. We identify the transition scale to homogeneity cor- 
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responding to which the value of 1 - lies within 10 ^ 

(Pandey 2013). This criteria is satisfied at 140/;^* Mpc for the 
sample 4. We also test if this result holds for galaxy sam¬ 
ples having sizes bigger than sample 4 and therefore probing 
larger length scales. We construct another two volume lim¬ 
ited samples namely sample 5 and sample 6 (see Table 4) 
which are larger than sample 4. We could not prepare the 
mock samples for sample 5 and sample 6 from the Millen¬ 
nium simulation due to the limitations arising out of the 
size of the simulation volume. However we analyze sample 5 
and sample 6 and find that regardless of their volumes and 
number densities the distributions become homogeneous at 
~ 140/i“* Mpc. This can be clearly seen in Figure 4 where we 
find that the value of the relative Shannon entropy of sample 
4, sample 5 and sample 6 agree within 0.02% at 140/i“^ Mpc 
conforming to homogeneity. Importantly there are ~ 400 cen¬ 
ters available at 140/i“* Mpc for both sample 5 and sample 
6 which reduces to 100 at 152/!“* Mpc and 166/!“* Mpc re¬ 
spectively. The observed offset between the curves below the 
scale of 140/i“* Mpc may indicate the presence of a sample 
size dependent finite size effect. To test this we applied the 
same cut in the upper redshift limit of sample 5 and sample 6 
as applied for sample 4 (Table 4). This ensures same volume 
and geometry for these samples as they have the same lower 
redshift limit and sky coverage. We reanalyzed the resulting 
samples and find that the offset between the curves per¬ 
sists with very little changes in the differences and all three 
curves finally merges at ~ 130/!“* Mpc (Figure 6). The offset 
between the curves primarily originate from the differences 
in the number density and clustering in the galaxy distri¬ 
butions. sample 5 and sample 6 have brighter magnitude 
limits and lower number densities than sample 4. Poisson 
noise induced inhomogeneities resulting from different num¬ 
ber densities plays an important role on small scales and 
become less important as the radius increases (Figure 2). 
Differences in inhomogeneities due to the differences in clus¬ 
tering may persist upto a large length scale. The fact that all 
three curves from sample 4, sample 5 and sample 6 reconcile 
at 140/!“* Mpc indicates that the differences in their inhomo¬ 
geneities cease to exist at this length scale. The length scale 
where all three curves merge reduces to ^ 130/i“* Mpc in this 
case (Figure 6). This results from the decrease in the avail¬ 
able number of centres at all length scales when the sample 
5 and sample 6 are trimmed. It may be noted in Figure 5 
and Figure 7 that the number of available centers reduces 
from ~ 1000 to 100 at ~ 120/!“* Mpc and ~ 130/!“* Mpc for 
sample 5 and sample 6 respectively when they are trimmed 
to match with the size of sample 4. 

It may be interesting to note here that an analysis of fil- 
amentarity in the Luminous Red Galaxy (LRG) distribution 
found that the filaments are statistically significant upto a 
length scale of 120/i“*Mpc (Pandey et al. 2011). Some re¬ 
cent studies (Einasto et al. 2015a, b) of the distribution of 
nearby rich galaxy clusters and superclusters in the SDSS 
find the evidence of shell like structures on 120/i“*Mpc. In 
different panels of Figure 3 we see that the mock galaxy sam¬ 
ples from the Millennium simulation also does not become 
homogeneous at least up to a length scale of 80/i“* Mpc and 
become homogeneous at ~ 100/!“* Mpc. We note that all the 
SDSS samples exhibit a greater degree of inhomogeneity as 
compared to their mock counterparts from the Millennium 


simulation. Nevertheless we can not deny that the complex 
issues of overlapping and confinement bias cast some uncer¬ 
tainty in the interpretations which can be resolved with yet 
larger samples with a reasonable number density. 

In future we plan to carry out analysis with the Lu¬ 
minous Red galaxy distributions (LRG) (Eisenstein et al. 
2001) and the Baryon Oscillation Spectroscopic Survey 
(BOSS) (Dawson et al. 2013) from SDSS which cover enor¬ 
mous volumes providing opportunities to probe much larger 
scales. Noticeably this would also enable us to perform 
the analysis with non-overlapping independent spheres even 
at sufficiently larger length scales. Analysis with non¬ 
overlapping spheres would have a remarkable advantage over 
the overlapping spheres as in the later the inhomogeneities 
are suppressed at each length scales making it difficult to 
detect them particularly at larger scales. This would help us 
to discern any inhomogeneities present even at the largest 
scales with nearly equal ease and settle the issue of cosmic 
homogeneity with an unprecedented confidence. 


5 ACKNOWLEDGEMENT 

Funding for SDSS-III has been provided by the Alfred 
P. Sloan Foundation, the Participating Institutions, the 
National Science Foundation, and the U.S. Department 
of Energy Office of Science. The SDSS-III web site is 
http://www.sdss3.org/. 

SDSS-III is managed by the Astrophysical Research 
Consortium for the Participating Institutions of the SDSS- 
III Collaboration including the University of Arizona, the 
Brazilian Participation Group, Brookhaven National Lab¬ 
oratory, Carnegie Mellon University, University of Florida, 
the French Participation Group, the German Participation 
Group, Harvard University, the Institute de Astrofisica de 
Canarias, the Michigan State/Notre Dame/JINA Partici¬ 
pation Group, Johns Hopkins University, Lawrence Berke¬ 
ley National Laboratory, Max Planck Institute for Astro¬ 
physics, Max Planck Institute for Extraterrestrial Physics, 
New Mexico State University, New York University, Ohio 
State University, Pennsylvania State University, University 
of Portsmouth, Princeton University, the Spanish Partic¬ 
ipation Group, University of Tokyo, University of Utah, 
Vanderbilt University, University of Virginia, University of 
Washington, and Yale University. 

REFERENCES 

Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2014, ApJS, 
211, 17 

Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, 
arXiv:1501.00963, Accepted to ApJS 
Amendola, L., & Palladino, E. 1999, ApJ Letters, 514, LI 
Baugh, C.M., Cole, S., Frenk, C.S. & Lacey, C.G. 1998, ApJ, 498, 
504 

Benson, A.J., Lacey, C.G., Baugh, C.M., Cole, S. & Frenk, C.S. 
2002, MNRAS, 333, 156 

Biswas, T., Notari, A., & Valkenburg, W. 2010, JCAP, 11, 30 
Bharadwaj, S., Gupta, A. K., & Seshadri, T. R. 1999, A&A, 351, 
405 

Bolejko, K., Celerier,M.-N., & Krasihski, A. 2011, Glassical and 
Quantum Gravity, 28, 164002 
Borgani, S. 1995, Physics Reports, 251, 1 


MNRAS 000, 1-12 (2015) 



12 Pandey, B. and Sarkar, S. 


Blake, C., Wall, J. 2002, Nature, 416, 150 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Buchert, T., & Ehlers, J. 1997, A&A, 320, 1 
Buchert, T. 2001, General Relativity and Gravitation, 33, 1381 
Clarkson, C. 2012, Comptes Rendus Physique, 13, 682 
Clifton, T., Ferreira, P. G., Land, K. 2008, Physical Review 
Letters, 101, 131302 

Cole, S., Aragon-Salamanaca, A., Frenk, C.S., Navarro, J.F., Zepf, 
S.F. 1994, MNRAS, 271, 781 

Cole, S., Lacey, C.G., Baugh, C.M.& Frenk, C.S. 2000, MNRAS, 
319, 168 

Coleman, P. H., Pietronero, L. 1992, Physics Reports, 213, 311 
Colles, M. et al.(for 2dFGRS team) 2001,MNRAS,328,1039 
Croton, D. J., Springel, V., White, S. D. M., et ah 2006, MNRAS, 
365, 11 

Dawson, K. S., Schlegel, D. J., Ahn, C. P., et ah 2013, AJ, 145, 

10 

De Lucia, G., Blaizot, J. 2007, MNRAS, 375, 2 
Fisenstein, D. J., Annis, J., Gunn, J. F., et ah 2001, AJ, 122, 
2267 

Fllis, G. F. R. 2011, Classical and Quantum Gravity, 28, 164001 
Finasto, M., Heinamaki, P., Liivamagi, L. J., et ah 2015, 
arXiv: 1506.05295 

Finasto, M., Gramann, M., Saar, F., et ah 2015, A&A, 580, A69 
Fixsen, D. J., Cheng, F. S., Gales, J. M., et ah 1996, ApJ, 473, 
576 

Gunn, J. F., Siegmund, W. A., Mannery, F. J., et ah 2006, AJ, 
131, 2332 

Gunn, J. E., Carr, M., Rockosi, C., et ah 1998, AJ, 116, 3040 
Guo, Q., White, S., Boylan-Kolchin, M., et ah 2011, MNRAS, 
413, 101 

Guzzo, L. 1997, New Astronomy, 2, 517 

Hawkins, E., Maddox, S., Cole, S., et ah 2003, MNRAS, 346, 78 
Hazra, D. K., Shafieloo, A. 2015, arXiv:1506.03926 
Hogg, D. W., Fisenstein, D. J., Blanton, M. R., Bahcall, N. A., 
Brinkmann, J., Gunn, J. E., & Schneider, D. P. 2005, ApJ, 
624, 54 

Kauffmann, G., White, S.D.M. Sz Guiderdoni, B. 1993, MNRAS, 
264, 201 

Kauffmann, G. Sc White, S.D.M. 1993, MNRAS, 261, 921 
Kauffmann, G. 1996, MNRAS, 281, 487 

Kolb, E. W., Matarrese, S., Sc Riotto, A. 2006, New Journal of 
Physics, 8, 322 

Kolb, E. W., Marra, V., Sc Matarrese, S. 2010, General Relativity 
and Gravitation, 42, 1399 

Kurokawa, T., Morikawa, M., Sc Mouri, H. 2001, A&A, 370, 358 
Larena, J., Alimi, J.-M., Buchert, T., Kunz, M., Sc Corasaniti, 
P.-S. 2009, Physical Review D, 79, 083011 
Martinez, V. J., Sc Jones, B. J. T. 1990, MNRAS, 242, 517 
Martinez, V. J., Sc Coles, P. 1994, ApJ, 437, 550 
Martinez, V. J., Pons-Borderia, M.-J., Moyeed, R. A., Sc Graham, 
M. J. 1998, MNRAS, 298, 1212 
Pan, J., & Coles, P. 2000, MNRAS, 318, L51 
Pandey, B. 2013, MNRAS, 430, 3376 

Pandey, B., Kulkarni, G., Bharadwaj, S., Sc Souradeep, T. 2011, 
MNRAS, 411, 332 

Paranjape, A. 2009, arXiv:0906.3165 
Penzias, A. A., Sc Wilson, R. W. 1965, ApJ, 142, 419 
Peebles, P. J. E. 1993, Principles of Physical Cosmology. Prince¬ 
ton, N.J., Princeton University Press, 1993 
Peebles, P. J. E. 1980, The large scale structure of the Uni¬ 
verse. Princeton, N.J., Princeton University Press, 1980, 435 
P-> 

Sarkar, P., Yadav, J., Pandey, B., Sc Bharadwaj, S. 2009, MNRAS, 
399, L128 

Scharf, C. A., Jahoda, K., Treyer, M., et ah 2000, ApJ, 544, 49 

Schwarz, D. J. 2002, arXiv:astro-ph/0209584 

Scrimgeour, M. L, Davis, T., Blake, C., et ah 2012, MNRAS, 3412 


Shannon, C. E. 1948, Bell System Technical Journal, 27, 379-423, 
623-656 

Smoot, G. F., Bennett, C. L., Kogut, A., et ah 1992, ApJ Letters, 
396, LI 

Somerville, R.S. Sc Primack, J.R. 1999, MNRAS, 310, 1087 
Springel, V., White, S. D. M., Jenkins, A., et ah 2005, Nature, 
435, 629 

Strauss, M. A., Weinberg, D. H., Lupton, R. H., et ah 2002, AJ, 
124, 1810 

Straumann, N. 1974, Helvetica Physica Acta, 47, 379 
Sylos Labini, F., Vasilyev, N. L., Sc Baryshev, Y. V. 2007, AScA, 
465, 23 

Sylos Labini, F., Vasilyev, N. L., Sc Baryshev, Y. V. 2009a, Eu¬ 
rophysics Letters, 85, 29002 

Sylos Labini, F., Vasilyev, N. L., Pietronero, L., Sc Baryshev, Y. V. 

2009b, Europhysics Letters, 86, 49001 
Sylos Labini, F., Sc Baryshev, Y. V. 2010, JCAP, 6, 021 
Sylos Labini, F. 2011a, Classical and Quantum Gravity, 28, 
164003 

Sylos Labini, F. 2011b, Europhysics Letters, 96, 59001 
Tomita, K. 2001, MNRAS, 326, 287 

Valkenburg, W., Marra, V., Sc Clarkson, C. 2014, MNRAS, 438, 
L6 

White, S. D. M. 1979, MNRAS, 186, 145 
White, S. D. M. Sc Frenk, C. S. 1991, ApJ, 379, 52 
Wilson, R. W., Sc Penzias, A. A. 1967, Science, 156, 1100 
Wu, K. K. S., Lahav, O., Sc Rees, M. J. 1999, Nature, 397, 225 
Yadav, J., Bharadwaj, S., Pandey, B., Sc Seshadri, T. R. 2005, 
MNRAS, 364, 601 

Yadav, J. K., Bagla, J. S., &; Khandai, N. 2010, MNRAS, 405, 
2009 

York, D. G., et ah 2000, AJ, 120, 1579 

Zehavi, L, Zheng, Z., Weinberg, D. H., et ah 2005, ApJ, 630, 1 
Zibin, J. P., Moss, A., Sc Scott, D. 2008, Physical Review Letters, 
101, 251303 

This paper has been typeset from a T^]X/IXI^jX file prepared by 
the author. 


MNRAS 000, 1-12 (2015) 


